Scale dependence of the quark masses and mixings: leading order 
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We consider the Renormalization Group Equations (RGE) for the couplings of the Standard Model 
and its extensions. Using the hierarchy of the quark masses and of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix our argument is that a consistent approximation for the RGE should be based on the 
parameter A = \Vus\ ~ 0.22. We consider the RGE in the approximation where we neglect all the 
relative terms of the order ~ A 4 and higher. Within this approximation we find the exact solution of 
the evolution equations of the quark Yukawa couplings and of the vacuum expectation value of the 

£-f) , Higgs field. Then we derive the evolution of the observables: quark masses, CKM matrix, Jarlskog 

invariant, Wolfenstein parameters of the CKM matrix and the unitarity triangle. We show that the 
angles of the unitarity triangle remain constant. This property may restrict the possibility of new 

£S) , symmetries or textures at the grand unification scale. 
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I. INTRODUCTION 

> . 

The solutions of the Renormalization Group Equations (RGE) determine the dependence of the physical parameters 
of the theory on the renormalization point [ESIESIEISHlElS IH1 EH • The set of these parameters for a given 
energy defines the theory that is equivalent to the one with the initial parameters. This property was used in the 
Standard Model in the search of the theory that is equivalent to the low energy Standard Model. This led to the 
hypothesis of grand unification when it was observed that 3 gauge coupling constants of the Standard Model converge 
to one value at the energy ~ 10 14 GeV, and to the reduction of the number of the parameters of the Standard Model. 
Further reduction of the number of the parameters is possible in the sector of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix |T^ . Il3l ] by looking for the textures or new symmetries. The CKM matrix appears in the Standard Model as 
a result of the transition from the quark gauge eigenstates to the quark mass eigenstates upon the diagonalization 
CIh, of the quark mass matrices. The quark mass matrices appear after the spontaneous symmetry breaking from the 
quark-Higgs Yukawa couplings. For this reason we consider the RGE for the quark-Higgs Yukawa couplings from 
which we obtain the evolution of the quark masses and the CKM matrix. 

In a recent paper 01 we began the analysis of the solutions of the one and two-loop evolution equations (obtained in 
a general quantum field theory) for the coupling constants, the quark Yukawa couplings (QYC) and the CKM matrix. 
In Ref. hj we systematically investigated the influence of the hierarchical structure of the QYC on the evolution 
of the CKM matrix by constructing the exact solution of the one loop RGE compatible with the observed hierarchy 
The most important result that we derived is that the CKM matrix evolution depends only on one universal function 
of energy which is a suitable integral that depends on the model (Standard Model (SM), Minimal Supersymmetric 
Standard Model (MSSM) and Double Higgs Model (DHM)). We showed that the evolution of the ratios of the 
eigenvalues (masses) of the up and down QYC, Y UtC /Y t and Y^/Yb depend on the same universal function as the 
CKM matrix. The eigenvalues of the QYC, Y u and Y c depend linearly on the corresponding initial values and their 
ratio Y u /Y c is constant while the functional dependence for Y t is non linear. The other remarkable result is that the 
diagonalizing matrices of the up quark Yukawa couplings are energy independent in the leading order. This means 
that the transformation (tpu)L,R — * (U u )l,r{iPu)l,r, will diagonalize the matrix of the up quark Yukawa couplings 
and it will stay diagonal upon the renormalization group evolution, and the evolution of the CKM matrix will be 
determined only from the evolution of the down quarks Yukawa couplings. This fact may simplify the model building 
based on the symmetries of the quark Yukawa couplings. 
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To achieve our present goal to obtain the precise results for the running of the parameters of the standard model 
and its extensions we will use the hierarchy and approximation scheme based on the parameter A = |V^ S | ~ 0.22 and 
establish accordingly the consistent procedure for the solution of the RGEs. We show that a one loop approximation 
is equivalent to neglecting in RGEs of the terms of the relative order A 4 and higher. This means that the terms of 
this type should also be neglected, if present, in the one loop RGEs. The inclusion of the terms of the relative order 
A 4 and A 5 requires the two loop part of the RGEs and for the terms of the order A 8 one has to include the three loop 
terms. 

Our aim is to systematically investigate the running of the observables of the currently available models up to the 
order A 5 which we plan to do in two steps. In the first step we consider the RGEs up to the order A 3 and use the 
fact that these equations can be analytically solved. In the second step we apply the modified perturbation calculus 
to obtain the evolution of the parameters with the precision A 5 . The methods and scope of each step are very wide 
apart so we will divide our analysis into two papers. In the present paper we establish the framework of our analysis, 
present exact solutions of the RGEs for the Yukawa couplings, vacuum expectation values and other observables and 
also show a graphical representation of the analytical results. In the forthcoming paper, using these results, we will 
develop the modified perturbation calculus which will be the basis of the derivation of the corrections up to the order 
A 5 to the analytical results from the first paper. 

The organization of the paper is the following. The Section [H] is devoted to the discussion of the hierarchy of 
the couplings and the observables in the RGEs and the introduction of the approximation scheme that we apply 
to the solution. In Section IIIII we discuss RGEs for the Yukawa couplings and the vacuum expectation values in 
the A 3 order and present the analytical solutions. The behavior of the universal energy function h(t) and of all the 
other functions that participate in the evolution of the physical parameters of the SM and its extensions is presented 
graphically from the mass scale mt to the Planck mass. In Section llVl with the exact solutions of the one loop RGE 
equations previously obtained, we solve, considering a new approach, the equations for the other observables (all the 
CKM matrix elements, the masses, unitarity triangle, the Jarlskog invariant J and the Wolfentstein parameters A, 
A, p, 77), we present their explicit analytical evolution which were not published before. Subsequently we represent 
graphically the RGE flow of these observables. In Section we draw the conclusions. Our new approach confirms 
the earlier results and also enables us to obtain the new ones, such as the differential equations for the diagonalizing 
matrices of the down sector in terms of which the CKM matrix elements are obtained. 



II. APPROXIMATION SCHEME 



In the approximate calculations the overall consistency is very important. This means that the terms that are 
smaller than those neglected should not be kept. The order of magnitude (expressed in the powers of A) of the 
components of the RGEs is 

g\ w 0.2 — A, <^«0.42~2A, g\ « 1.5 ~ 1, y u yt ~ 1, 

/ \ 2 1 (1) 

VdVa~ — -1-13A 5 , — -r, = 6.3 • 10- 3 « 2.7A 4 , Xh ~ 0.44 ~ 2A. 
\m t J (4tt) 2 

Here gi are the gauge coupling constants, y u and yd are the matrices of the Yukawa couplings and Xh is the Higgs 
quartic coupling. 

The mathematical structure of the RGEs is the following 

f = (4^ {"* ^ + W?^ {glVuVl VdVl) + " ' } 91 (2a) 
= J^y tedCs*. VvVl, Vdvl) + j^Pu,d(9l, VuVl, Vdv\, Ah) + • • ■ I y u ,d (2b) 

~dl = (i^) 1 VdV ^ + (i^j 1 ^ VdV ^ + " ' } V ^ 



dX 
~dt 



In Eqs. J2J) t = ln(E/fi) is the energy scale parameter, /1 is the renormalization point, v is the the vacuum expectation 
value of the Higgs field, (3^ denote the n-loop (3 functions for the RGEs which are homogenous polynomials of the 
indicated variables and the dots indicate the omitted three and more loops contribution. The approximation scheme 
for the RGEs consists in the approximation of the terms inside the braces on the right hand side of the equations 
keeping only the terms of a given relative order. 
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From Eqs. J5J and (JTJ one can see that the leading terms of the one loop contribution are of the order 1, because 

P {1) {gl...)-i (3) 

and for the two loop contribution 

' ^(^...)~7^~^. (4) 



{At:) 2 ' v ™' ' (4tt) 

From Eqs. © and (@J we thus see that in the one loop approximation we should neglect all the terms of the order 
A 4 and higher. From the hierarchy given in Eq. Q it follows that the one loop approximation is equivalent to the 
neglecting the two loop term and also putting yd,y\ = on the right hand side of relatively to the terms of the 
lower order. Such equations, as it turns out, can be analytically solved. 

The next order term is of the order A 4 . Such a term is present only in the two loop contribution and it does not 
introduce important qualitative modification to the running of the observables and it will be discussed together with 
Eq. I|2djl and the A 5 contribution in the forthcoming paper. RGEs for the order A 4 and higher form a system of coupled 
nonlinear equations which are difficult to solve analytically, so for their analysis one has to use other methods. 

III. ONE LOOP EQUATIONS AND SOLUTIONS FOR THE GAUGE COUPLING PARAMETERS, THE 
YUKAWA COUPLINGS AND THE HIGGS VACUUM EXPECTATION VALUE 

In this section we discuss the one loop RGEs for the gauge, Yukawa couplings and the vacuum expectation values. 
This means that the precision of all the expressions is up to the order A 3 . We indicate this fact in the final form of 
the equation for each observable. The structure of the one-loop RGE for the gauge coupling parameters gk is 



dt (4tt) 2 



hgi + 0{\*). (5) 



Here the coefficients (h,b 2 ,b 3 ) are equal (41/10,-19/6,-7), (21/5,-3,-7) and (33/5,1,-3) for the SM, DHM and 
MSSM, respectively. The solution to this equation is derived directly 



The one-loop RGE for the Yukawa couplings y Ut d are 

dy u ,d 



dt L(4tt) 



1 -d) 



Pu,d 



yu.d + 0(\ A ). (7) 



The /3$pj has a hierarchical structure based on the parameter A and has the following form 1 (with lepton dependence 
suppressed) 

PP = a\(t) + aUy u yt) + aiTr(y u yt) + a\(y d yl) + 4^{y d y% i = u,d. (8) 

Its structure allows the systematical solution of the RGE with A as the expansion coefficient. The approximate form 
of the equations for the quark Yukawa couplings, neglecting all the terms of A and higher is the following 

dVu 1 r..«™ , „t , „¥T>r„ „t 



(U [4 ^ L aUt) + a^y u yt + a^Tr(y u yt)]y u + 0(X% (9) 
1? = j^[*t(t)+aiy u yi + aiTr(y u yt)}y d + 0(\% (10) 



1 Notice that the form of the function is compatible with Eq. 12bt . i.e. it is the homogenous polynomial of order 1 of the indicated 
variables. 
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where 

«i(t) = -(ci.9? + c 2 g 2 2 + c 3 53 2 ), «2=y, "3=3, 
<*i(t) - -(^9l + c' 2 9l + c' 3 gl), a d 2 = y, ag = 3a, 

(a,6,c) are equal to (1, 1, -1),(0, 1, 1/3), (0,2,2/3); (ci,c 2) c 3 ) are equal to (17/20,9/4,8), (17/20,9/4,8), 
(13/15,3,16/3) and (ei,4,c(j) are equal to (1/4,9/4,8), (1/4,9/4,8), (7/15,3,16/3) in the SM, DHM and MSSM, 
respectively. 

The transformation from the quark gauge states to the physical states requires the diagonalization of the Yukawa 
coupling matrices y u and y d with the biunitary transformations 

Diag(y u ( d ),F c(s ),F 4(b) ) = {U u ( d )) L y u ( d) {U u[d) % (11) 

where (U Ut d)L,R are the corresponding unitary diagonalizing matrices and Y u ^, Y" c ( s ) and Y t ^ are the eigenvalues of 
the Yukawa coupling matrices y u r d ) |15|| . As it is well known, the diagonalizing matrices generate the flavor mixing in 
the charged current described by the CKM matrix 

V = (U u ) L (U d )l 

The same biunitary diagonalization transformation also permits the exact solution of Eq. (0 which arises from the 
fact that the diagonalizing matrices (U u )l.r for the one loop Eq. © do not depend on energy ^4[. It then follows 
that y u (t) from Eq. (|5J has the following representation 

Vu{t) = (U u )[Y u (t)(U u ) R , Y u {t) = Diag(Y„, Y c ,Y t ) (12) 

The whole dependence on t in Eq. (|12|) is contained only in the diagonal matrix Y u (t). The diagonal elements of 
Y u (t) satisfy the system of differential equations that follows from Eq. J^J 



dY 

11 1 u,c 

dt 
dY 



(47T) 



: {aUt) + a u 3 Y t 2 }Y U:C + 0(X A ), 



~JT = 77TT7{«rW + K + «3)^ 2 }^+0(A 4 ) 



dt 



(4tt) 



(13a) 
(13b) 



Eq. (|13b|l decouples from Eqs. (|13all and can be solved independently of the other equations. When the solution for 
Y t (t) is known then Eqs. H13a|) become linear and can also be solved and eventually the solution of Eqs. (|13|) is 



Y u , c (t) = Y u , c {to)VHt)h%? (t) , 
Y t (t) = Y t (t )Vmhif 2+a *\t). 



where the functions r(t) and h m {t) are 

2 



r(f) 



h m (t) 



exp 



exp 



_(4tt) 2 



a%(r)dT 



k=3 

n 

fc=i 



9l (to) 
9l(t) 



(14a) 
(14b) 



(15) 



(16) 



It is worth mentioning that the function h(t) = [h m (t)]" 2 appears in the evolution of the V matrix and moreover this 
is the only dependence on t of the V matrix. 

The next step is the determination, from Eq. I|10|) . of the running of the down quark Yukawa couplings. The 
substitution 



y d (t) = (U u )[W(t) 



(17) 
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transforms Eq. I|10[l into the following equation 

= jhf + a " {YU{t))2 + a " Tr (yu (<))2 } w ' 

with the matrix jaf(i) + a$ (V u (i)) 2 + "fTr (V"(i)) 2 } being diagonal. 
This allows to solve Eq. (|18[) explicitly and the solution for W(t) reads 

W(i) = (r'(i)) 1 ' 2 }d{t) ■Z(t)-W(t ) 

where Z(t) is the diagonal matrix 

Z(t) = Diag(l, l,h(t)). 

Now, with the help of Eq. (|17|l one gets the explicit one loop running of the Yukawa couplings yd(t) 

y d (t) = ^/?{t)h^(t)(U u )[Z(t)(U u ) L y a d , y° d = y d (t Q ), 

where 

fc=3 



(18) 



r'{t) — exp 



n 



k=l 



9k fa) 



The approximated one-loop RGE for the Higgs vacuum expectation value (VEV) v u , d is the following 

dv ud 1 



dt 



(4tt) : 



(19) 
(20) 

(21) 
(22) 

(23) 



where v u is the VEV for the up quarks and v d for the down quarks 2 . Note that for the SM v — v u — v d - The 
functions a v i u,d (t) and the coefficients are equal 



a 3 ° = -3, a," = -3, a,' = 0, 



(24) 
(25) 



and the constants (c'[ {u d) , 4 {u d} ) are equal to (-9/20,-9/4), (-9/20,-9/4) and (-3/20,-3/4) for the SM, DHM 
and MSSM, respectively. 

To solve Eq. 1123(1 for the VEV we divide both sides of the equation by the corresponding VEV and obtain after the 
integration 



1 



dv u ,d 

t Vu,d (4tt) 2 J t() 



^(t) + al^Tr(y u yi))dt. 



(26) 



The left hand side of Eq. (|26j) is equal to In [v u ,d (t) /v u ,d {to)} and on the right hand side we have the sum of two 
integrals. The first integral can be explicitly integrated using Eq. © and we introduce the function r" u (t) 



r " u M = ex P 



(4tt) 2 



a v ^ d (T)dT 



n 



fffc(*0) 



c fc(n,d) 



and the second integral is calculated using Eq. i|16|) 



(4tt) 2 



(4tt) 2 



Y t 2 (r)dT = al ud \Yih m (t) 



(27) 



(28) 



Co 



We consider only the case where tan/3 = v^/v u ~ 1. Large values of tan /3 require other treatment and will be considered elsewhere. 
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The final form of the VEV is thus equal to 

v u ,d (*) = v u>d {t ) y/r» u Jt)h% U ' d (t), < d = v u4 (t ) . (29) 

The solutions of the renormalization group equations for the gauge coupling constants, Eq. JJjJl, for the quark Yukawa 
couplings Eqs. 1)140. H21|) and l)43|) and for the vacuum expectation values Eq. 129fl form the complete set of the 
evolution functions from which one can obtain the renormalization group flow of all observables related to quarks: up 
and down quark masses and the CKM matrix. In the next sections we will analyze these observables but let us notice 
here that the evolution is described by the following functions of energy 

r(t),r'(t),r'ljt),h m (t),h(t). (30) 

The dependence of the observables on these functions will be discussed in the next sections. Here in Figs. ^ and [5] we 
show the functional dependence of the functions in Eq. 1)30(1 for three models (SM, DHM and MSSM) to be able to 
see what is their influence on the observables and how they depend on the model. Notice that in all the figures we 
choose as the renormalization point the mass of the top quark nit = 174.3 GeV 16]. In such a way the functions in 
Eq. 1)30)1 and observables are independent of the quark mass thresholds. The extensive discussion of the thresholds 
effects is given in Ref. Q. 



IV. EVOLUTION OF THE OBSERVABLES 



The solution of the renormalization group equations presented in Section II I II allows the analysis of the evolution 
of all the observables related with the Yukawa couplings and the Higgs vacuum expectation values, i.e. the quark 
masses and the Cabibbo-Kobayashi-Maskawa matrix. The results of this section have been obtained from the explicit 
solutions from the previous section. Their validity and precision are therefore the same, i.e. the terms of the order A 4 
and higher have been neglected. We start with the analysis of the quark masses presenting first the analytical results 
and then showing the corresponding graphs. The same type of the analysis will be also applied to the CKM matrix. 



A. Quark masses 



The quark masses after the spontaneous symmetry breaking are equal to 

tm = ^=Y t (31) 

where Yi are the eigenvalues of the corresponding Yukawa couplings and Vi is the vacuum expectation value of the 
Higgs field. For the theories with one Higgs doublet (SM) there is one Higgs vacuum expectation value and for two 
Higgs doublets (DHM, MSSM) there is one VEV for the up quarks and another for the down quarks. 



1. Up quark masses 



The evolution of the eigenvalues (Y u , Y c , Y t ) for the up quarks is given in Eqs. I|14|) and the evolution of the VEV's 
is given in Eq. 1)29( 1 . Using Eq. (|31|l we thus obtain 

< (*) = = V ^Y?{t )^r{t)r»^)h% (t) = m^{t )^ r{t)r>> u {t)h% (t) , (32) 

and the power Kf is equal 

K u,c = "3 + a 3 7 K t = a 2 + a 3 + a 3 ■ I 33 ) 

In Figs. |3 (a) and 0(b) we show the running of the ratios of the up quark masses mi(t) / 'm^o) i = u,c,t for three 
models SM, DHM and MSSM. 
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1.10 



1.05 



1.00 



0.85 



0.80 




t = In ( E/m t ) 



FIG. 2: The scale dependence of the ratio h(t)/h(mt) 



2. Down quark masses 

The evolution of the eigenvalues of the down quark Yukawa couplings is more complicated because the diagonalizing 
matrices are also running. If we write the matrix y d in the form 



y d = (U d ) ' Y d (U d ) , 



then Eq. I) Kill becomes 



dt {At:) 2 
which after some simple manipulations becomes 
rd\t 



U d ^^Y d + + Y d ^ (U d R y = (Ul) {af(t) + a d yj u + afErfo.yt)} y 



(4tt) : 



(34) 



(35) 



(36) 



(37) 



where the matrix M vv is introduced for notational simplicity. 

The matrices U d ^$^- and ^f- (U d ) 1 
are purely imaginary. The diagonal elements of the matrix Yd are purely real and the matrix M vv is hermitian, so 
after taking the real part of the diagonal elements of Eq. (|36|l we obtain 



L — — an d d t (Ur) are antihermitian (this follows from unitarity) so their diagonal elements 



( dY d 



1 



(M vv Y d ) 



(38) 



V dt J u (4tt)2 
The matrix [M vv )^ has the simple form in our approximation. 

{M m ) l3 = (v f af(t)V + V*c4Y*V + yt^Tr(F u 2 )y) . = {af(t) + a d Y?) % + a d Y?V^V tJ + 0(A 4 ). (39) 
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The last step in derivation of l|39|) follows from the approximation based on hierarchy 



(Y u ) 



X 35 Y,. 



Finally we obtain the following equations for the eigenvalues of the down quark Yukawa couplings 

rlV d 1 - 

W-T,',' ^ | J A^ r O -r% A-w rO. I •* rA /H /" \ 4 \ 



dt 



(At:) 2 









{<4 + ^y t 2 


Vu 


2 + 4Y t 2 ) 



Using the hierarchical properties of the V matrix 

2 



V, 



hi 



A 4 , 



V 



lb 



= 1, 



(40) 



(41) 



we obtain the following equations for the eigenvalues of the down quarks Yukawa couplings 



The solution of Eqs. l(4*2|) reads 



dY d , s 


1 


dt 


(4tt) 


dYb 


1 


dt 


(4tt) 




Yd At) 




Y b (t) 



(42) 



(43) 



where the function r'(t) is given in Eq. These results complement previous results shown in Sec. IIIII Now, 

Eqs. Q43[l together with the VEV evolution Eq. (|29|) gives the evolution of the down quark masses 



mf(t) = ^Yfit) = V J^Y l d {t )^r\t)r^ d {t)h^ (t) = m f{t )^ r> {t)r>' d {t)h« f (t) 



V2 

there the powers Kf are equal to 



V2 



Ki = ai + ai + al d . 



(44) 



(45) 



Eqs. I|44|l give the analytical form of the down quark mass evolution. The evolution of the ratios of the down quark 
masses for three models SM, DHM, and MSSM is shown in Fig. 01(c) and Fig. El (d). 



B. Cabibbo-Kobayashi-Maskawa matrix 



The other set of observables is related with the CKM matrix. These observables include the absolute values of the 
matrix elements of the CKM matrix, Jarlskog's parameter J ^tJi Wolfenstein parameters A, A, p, r\ 0,^1 an d the 
angles of the unitarity triangle |2lt |22| . 

The evolution of the CKM matrix is simpler than that for masses because the CKM matrix depends only on the 
left diagonalizing matrices of the biunitary transformations for the Yukawa couplings and as was shown in Ref. [l4T | 
the evolution of CKM matrix depends only on one function of the energy h(t). This fact implies that there exist 
correlations between the evolution of the various elements of the CKM matrix. 

In the following we shall first find the evolution equations of the absolute values of the CKM matrix elements and 
afterwards analyze the J parameter and the unitarity triangle. 



1. Evolution of the absolute values of the CKM matrix V. 

Our starting point is Eq. (|36|l . We take now the imaginary part of the diagonal elements and the full off-diagonal 
elements of Eq. I|36|l , and using the form of the matrix M vv in Eq. I|39|) we obtain 



dt 







(46) 
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and 



(47T) 



(47) 



The hermitian conjugate of Eq. I|47|l reads 



yf (ui) 



d(ui) r 



'd(U d R ) 



i + 3- 



(48) 



Now multiplying Eq. (|4TJl by and Eq. (|4*8)) by and summing these equations, we obtain after some simple 
manipulations the differential equation for the evolution of the right diagonalizing matrix of the down sector 



'd(U d R ) 



2Y d Yf 



{Uf} 



(49) 



Eq. I|49|l gives the relation between the evolution of the left and right diagonalizing matrices and it is the key relation 
that permits the derivation of the evolution of the left diagonalizing matrix of the down quark Yukawa couplings. 
Inserting Eq. (|49|l in Eq. I)47|) wc obtain the evolution equation for the U d matrix 



(yf) +m 



dt 



(4tt) 



T a d 2 Y t 2 V*V tJ 



(50) 



We will now convert Eq. I|5U|I into equations for the CKM matrix elements. To this end we first notice that within 
our approximation (neglecting all terms of the order A 4 and higher) 



-sign(i - j) + 0(A 4 ), 



(51) 



where sign(i — j) — =pl for i ^ j. 

Next we know that the up quark diagonalizing matrix does not depend on the energy. Using the definition of 



the CKM matrix V = U% [Uf) ] we obtain 



dt 



Now from Eqs. I|50(l . I|51|l and (|52(l the following equations are obtained 

1 



dt 



(4tt) : 



;sign(i-j)a^Y t 2 V t *V tj . 



(52) 



(53) 



Eq. I|53l) is valid only for i j but from the property that the matrix (v^dV/dtj is antiunitary and its diagonal 
elements are purely imaginary, we obtain equations for the absolute value of the elements of the CKM matrix 



V cd 



V cd 



dt 



Vtd 



dt. 



(47T) 



V td 



(47T) 



V td 



(54a) 



(54b) 



V ub 



V ub 
~dt 



d v 2 



(47T) 



l^t 



v tb 



(54c) 



V, 



v cb 



cb 



dt 



d v 2 



(47T) 



2°& Y t 



V tb 



(54d) 
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and from the unitarity of the CKM matrix we obtain the evolution of the remaining CKM matrix elements. 
To solve Eqs. I|54fl let us notice the following relation 

1 dh 2 



a*Y 2 = TTT ~ 



(4tt) 2 2 * h 2 dt ' 

where h (t) = [h 7n (t)] a2 and h m (t) is the function introduced in Eq. (|16fl . Using relation (|55fl we obtain from Eq. O 



(55) 



|Vi rf (t) 



From Eqs. (|54a|l and (|54b() we obtain 



IK/ 



/ l 2 + (1 _ /l 2 ) | Fi 0|2 



where V? = % (to) . 



\V cd (t) 



h 2 \V c ° d \ 2 



h 2 + {1 - h 2 )\V t ° d \ 2 

Using the relation 1 — \V tb \ 2 = \V ub \ 2 + \V cb \ 2 we derive from of Eqs. I54cf) and (|54d|) the following result 3 



\v tb (t)\ 2 = 



h 2 \V° b \ 2 



l + (h 2 -l)\V t l\ 



and this result yields with the help of Eqs. (|54c|l and (|54dl ) the evolution of the | V^b(t) | 2 and |V^b(i)| 2 : 

\V ub (t)\ 2 



V cb (t) 





v° h 


I 2 


l + (h 2 - 


m° b \ 2 






2 




K° 6 





i + (tf-i)|vsi a 

The remaining elements of the CKM matrix are obtained from the unitarity relation 



V ud (t) 



h 2 



ud 



h 2 + (l-h 2 ) v t ° d 



(56a) 



(56b) 



(56c) 



(56d) 



(56e) 



(56f) 



V us (t) 



V cs (t) 



Vts (t) 



v° 

us 



+ {h 2 1) 



2 ,2 
= V° 



K° 



h 2 - 1 



(h 2 i) 



[ IK?/ 


IK / 


IK?/ 


IK/ 


l + (h 2 - 


1)IK°/ 


h 2 + (1 - 


^ 2 )IK /. 


IK / 


K / 


IK / 


K/ 


_l + (h 2 - 


1)IK°/ 


/i 2 + (i - 


^ 2 )IK / 


[lK°/(i- 


- IK /) 


IK°/(i- 


- IK/) 1 


l + (h 2 - 


1)IK°/ 


h 2 + (1 - 


^ 2 )IK / 



(56g) 



(56h) 



(56i) 



Eqs. (|56|l represent the RG evolution of the CKM matrix. The absolute values of the CKM matrix elements do not 
depend on the parameterization of the CKM matrix. Notice that in agreement with the Theorem 1 in Ref. [l4| the 
evolution of the CKM matrix depends only on one function h(t) of energy. It is for this reason that we consider h(t) 
as an universal function of energy. 



3 Note that the equation for the evolution of |Vtb| 2 given in Ref. Q, Eq. (39d) has a missing h 2 in the numerator. 
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2. Jarlskog Invariant 



Jarlskog parameter J is defined as [l7j 



J = 3 



VudV cs V* s V* d = 9 V ud V t6 VM =9f(I>), 



(57) 



and it is nonvanishing if the CKM matrix is non CP invariant. 
Let us now consider the following expression 



dhxD dln [ v ndVtbV ub V td Id. 1 d- 1 d.* 1 d Tr , 

V ud dt Vtb dt v: b dt 



dt 



dt 



V td ^ 



(58) 



Now using Eq. (|53l) and the property that the matrix V' ^ is antihermitian (and has imaginary diagonal elements) 
we obtain 



dlnD 



dt 



= ~M Y t 







2 








( 


Vtb 




v td 


2 ) « -2«^ t 2 





and thus the Jarlskog invariant fulfills the equation 

din J 



dt 



which gives using Eq. H56c[l 



2a d 2 Y 2 



In — = ln 

Jo 



Vt 



Lb 



1 dh 2 ~ 2 



1 + {h 2_ lW l2 



or 



J 



Jo 



l+(fc»-l)|V$P 



Jo 



(59) 



(60) 



(61) 



(62) 



The evolution of the Jarlskog invariant is given in Fig. 0] 

3. The Wolfenstein Parameters 



Eqs. I|56|l can be written in an approximate form using the hierarchy of the CKM matrix, given by the Wolfenstein 
parameterization [I8I.IT9I] 



V 



l-\X 2 



-A 

AA 3 (l-p-^) -AX 2 



A AX 3 (p-irj) 
AX 2 
1 



1-iA 2 



P = P 1 



Now neglecting all the terms of the relative order A 4 and higher, we obtain 



Vao 



1 - 



v u ° d 



v cd 



v c ° d 



Vtd 



v C! 



1 Id 



v ut 



v° 



v t ° s 



v° 

us 

1-h 2 
h 2 



V ut 



r\ = 77 1 - 



VS. 



A 2 



h 2 



V 

• /■ 



v cl 



v°, 

CO 



1 



v tb 



T-0 

v lb 



h 2 



(63) 



(64) 



From Eqs. (|64|l immediately follows the simple evolution of the Wolfenstein parameters 

A (t) = , and A, p, rj are invariant. (65) 

Notice that the dependence of the CKM matrix and Wolfenstein parameters on the renormalization scheme is given 
in Ref. E3. 
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5 10 15 20 25 30 

t = In ( E/m ( ) 



FIG. 4: The scale dependence of the ratio J(t)/J(mt) for the Jarlskog invariant. 



4- The Unitarity triangle 

The unitarity triangle plll^ is obtained from the scalar product of the first column of the V matrix by the complex 
conjugate of the third 



v V dV* b + v cd v; b + v td v; b = o 

and then by rescaling it in such a way that the lengths of the sides of the resulting triangle are equal 



(66) 



*-<>-t>5 



V,,, 



cb 



*4 



iii 



v, 



cb 



1. 



(67) 



From Eqs. (|64|l and (|65J) it immediately follows that R b and Rt are invariant upon the evolution which implies that 
the unitarity triangle is also invariant. Thus the complex phases of the CKM matrix elements Vtd and V ub (angles (3 
and 7) are also invariant up to the order A 4 . 

V. CONCLUSIONS 

In this paper we analyzed the solutions of the RGE for the quark Yukawa couplings, for the Higgs VEV's and also 
the evolution of all the observables that follow from them. The results depend on the model dependent functions 
given in Eq. J20J- The most interesting is the universal function h(t) because only on this function depends the 
evolution of the CKM matrix, Eq. I|56|) or (|64|l . The running of the absolute values of the CKM matrix elements and 
the invariance of the unitarity triangle angles implies that the evolution of the CKM matrix is remarkably simple 



V°, V°, V\ 

v ua y ua y uo 

v? 



cd 



Iff, Q ^^/l 



cb 



t>0 v o v o 
v td v ts v tb 



( 



v u ° d f& 



v° 



v° h \ 

ub 
cb 

h 



(68) 



We observe (see Fig- El) that the function h(t) is decreasing for the SM and increasing for the DHM and MSSM. This 
results in a qualitatively different evolution of the matrix elements and the observables related to the CKM matrix for 
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the SM in comparison to the DHM and MSSM. An important result is the dependence on it of the Jarlskog invariant 
which is shown in Fig. 0] We thus see that the CP violation is enhanced with increasing energy in the SM while it 
decreases for the DHM and MSSM. 

The evolution of the unitarity triangle is also very important. We have shown that the angles of the unitarity 
triangle remain constant upon the evolution. The invariance of the angles of the unitarity triangle is very significant 
because it means that at the grand unification scale the angles are the same as at the low energy so if there is a 
symmetry at the grand unification scale then it has to predict the low energy angles of the unitarity triangle. This 
strongly constrains the possible symmetries or textures at the grand unification scale. 

As we discussed in Section H] our results are approximate and we kept the terms up to the order A 3 . In the next 
order one has to keep the powers up to the order A 4 . In this approximation the results are qualitatively the same 
and the only corrections are small modifications of the functions h m (t) and h(t). The next order, A 5 , is significantly 
different and it will be discussed elsewhere. 
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